Spin susceptibility in bilayered cuprates: resonant magnetic excitations 
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We study the momentum and frequency dependence of the dynamical spin susceptibility in the 
superconducting state of bilayer cuprate superconductors. We show that there exists a resonance 
mode in the odd as well as the even channel of the spin susceptibility, with the even mode being 
located at higher energies than the odd mode. We demonstrate that this energy splitting between the 
two modes arises not only from a difference in the interaction, but also from a difference in the free- 
fermion susceptibilities of the even and odd channels. Moreover, we show that the even resonance 
mode disperses downwards at deviations from Q = {■k,tt). In addition, we demonstrate that there 
exists a second branch of the even resonance, similar to the recently observed second branch (the 
(5*-mode) of the odd resonance. Finally, we identify the origin of the qualitatively different doping 
dependence of the even and odd resonance. Our results suggest further experimental test that may 
finally resolve the long-standing question regarding the origin of the resonance peak. 

PACS numbers: 71.10.Ca,74.20.Fg,74.25.Ha,74.72.-h 



I. INTRODUCTION 



Magnetic excitations in the high-temperature super- 
conductors are of fundamental interest. While it is cur- 
rently still a topic of intense debate whether a contin- 
uum of magnetic excitations is responsible for the oc- 
currence of superconductivity in the cuprates, the feed- 
back effect of d2;2_.j^2-wave superconductivity on the mag- 
netic excitation spectrum has been well established in the 
context of the "resonance peak". This peak has been 
observed by inelastic neutron scattering (INS) experi- 
ments in three different families of the high-temperature 
superconductorsii^. The doping dependence of the peak 
frequency, rij.es(Q), the downward dispersion of the res- 
onance, which tracks the momentum dependence of the 
particle-hole continuum, and the emergence of a second 
resonance branch further away from Q are all consistent 
with the idea that the resonance peak is a particle-hole 
bound state (i.e. a spin exciton) below the particle-hole 
continuum. According to theory"', this excitonic reso- 
nance is a fundamental property of a (1^2 _y2 supercon- 
ductor. (For a review of other theoretical scenarios, see 

Recent INS experiments in overdoped YBa2Cu306+i 
(YBCO) revealed the formation of two resonance modes 
that differ by their symmetry with respect to the ex- 
change of adjacent copper oxide layers^. The original 
resonance mode observed in the bilayer cuprate possesses 
an odd (o) symmetry while the new one exhibits an even 
(e) symmetry. The frequency of the even mode is larger 
while its intensity is smaller than that of the odd mode. 
Moreover, while the doping dependence of the odd mode 
is non-monotonic and roughly follows 0°^^ ~ 5k bT^, 
the frequency of the even mode increases monotonically 



with decreasing dopingiiii^. 

The splitting in energy between odd and even reso- 
nances has been analyzed theoretically in the past within 
the random phase approximation (RPA)^^, and has been 
attributed to the difference in the strength of the resid- 
ual interaction leading to the bound state. The larger the 
interaction, the more the resonance is shifted downwards 
from the lower edge of the particle-hole (ph) continuum. 
How such a difference in the interaction can easily be 
seen in the t — J model, where the interactions in the 
even and odd spin channels are given by 



(1) 



with J|| , > being the in-plane and out-of-plane ex- 
change interaction, respectively. Thus Jo > Jei and the 
odd resonance occurs at a lower energy than the even one. 
Moreover, since the even mode lies closer to the ph con- 
tinuum its intensity is lower than that of the odd one. 
These two theoretical resultai^i^^ii^ are in good agree- 
ment with the experimental observation o^i^i^^'^^ . 

In this article, we address three issues which have not 
yet been considered in earlier studies on the spin reso- 
nance in bilayer systems. First, we argue that the differ- 
ence between the even and odd modes comes from two 
factors. One is the difference in the interaction, which 
was taken into account in earlier studies, another is the 
difference in the free-fermion susceptibilities of the even 
and odd channels which so far has been neglected. We 
show that the two factors are generally comparable to 
each other and depend on the same combination of pa- 
rameters. Numerically, the difference in the interactions 
leads to a larger splitting between the even and odd reso- 
nances than the difference between the even and odd free- 
fermion susceptibilities. Second, we extend our previous 
analysis of the odd resonance's dispersior^i^ to the even 
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channel, and show that the even resonance mode also 
disperses downwards at deviations from Q. Moreover, 
we show that the downward dispersion of the even mode 
is more parabolic than that of the odd channel. Third, 
we demonstrate that there exists a second branch of the 
even resonance, similar to the recently observed second 
branch (the Q*-mode^^) of the odd resonancefi2ii^. We 
show, following the approach of Ref.^^, that in the even 
channel, this second branch is much narrower in energy 
than in the odd one. These results suggest further exper- 
imental test that may finally resolve the long-standing 
question regarding the origin of the resonance peak. 

Finally, we analyze the doping dependence of the even 
and odd resonances. In the overdoped region, both 
modes decrease due to a decreasing superconducting gap. 
In the opposite limit of zero doping even and odd res- 
onances very likely evolve into the acoustic and opti- 
cal spin wave modes of the bilayer Heisenberg antifer- 
romagnet. We show, however, that, while plausible, the 
crossover from one regime to the other cannot be ob- 
tained within a simple RPA scheme chiefly because of the 
incorrect doping dependence of the free-fermion suscep- 
tibilities: the real part of both, the even and odd suscep- 
tibility decreases with decreasing doping at half-fiUingi^. 
This behavior is a direct consequence of the fact that the 
even susceptibility diverges at the van-Hove singularity, 
and the odd susceptibility possesses a maximum near the 
van-Hove point. 

The rest of the paper is organized as follows. In Sec.HIl 
we introduce our theoretical model and discuss the origin 
of the splitting between the even and odd resonance at 
Q = (7r,7r). In Sec. |III]we present the dispersion of the 
two resonances away from Q and show that a Q*-mode 
also arises in the even channel. In Sec. II VI we discuss the 
doping dependence of the resonances. Finally, in Sec. |V] 
we summarize our results and conclusions. 



II. EVEN AND ODD RESONANCES AT 

Q = ('r,7r) 

The coupling between two Cu02-planes in a unit cell 
of YBCO is described by the interlayer hopping matrix 
element i±(k) — jt± {cosk^ — cos ky), Ref. l20!. This 
coupling leads to the formation of bonding (6) and an- 
tibonding (a) energy bands whose dispersion are given 
by 



in the two layers via 



LL. 



'k'' ~ — 2t (cos fc:r -I- cos fcj^) 



±— tj^ {coskx — COS ky)^ — fi 



At' COS kx COS ky 



(2) 



with t = 250mcV, t' /t = 0.4, tj./* = 0.2, and being the 
chemical potential (these parameters provide a good fit 
to the Fermi surface of the bilayered Bi2Sr2CaCu2 08+5 
(BSCCO)^). The resulting Fermi surfaces for the bond- 
ing and antibonding bands are shown in Fig. [TJ 

The bonding and antibonding creation and annihila- 
tion operators are related to the fermionic operators, ci^2 



_ Ci+ C2 Ci — C2 



V2 



72 ' 



(3) 



It is also convenient to introduce even and odd compo- 
nents of the spins at site i, which are given by 



Sl(i) + S2(i) 



\ (4,a(i)0-a,/3Ca,/3(i) +4„(i)cra,/3Cb,/3(i)) 



Si(i)-S2(i) 



\ (4,a(i)0-a,/3Cb,/3(i) +cL (i)cra,/3Ca,/3(i) 



Se(i) 



So(i) 



(4) 

The experimentally measured susceptibility is related to 
the even and odd susceptibilities, =< SeSe > and 
X" =< SoSo > vie^ 

X(q,c.) = x^(q,c.)cos2^-fx°(q,^)sin2^, (5) 

where d is the separation between the layers. For non- 
interacting electrons, the susceptibilities in the even and 
odd channels are given by Xo = Xo" + Xo" and Xo = 
Xo^ + Xo"j respectively, where Xo° ^^'^ Xo' represent in- 
traband particle-hole excitations, and Xo' and Xo° rep- 
resent interband excitations. The free-fermion suscepti- 
bilities in the superconducting state at T = are given 
byi^ 




FIG. 1: (color online) Calculated Fermi surface(FS) for 
the bilayered cuprates as obtained from Eq.(12]l. The ar- 
rows indicates the transition between bonding-bonding (66), 
antibonding-antibonding (aa) and antibonding-bonding {ab, 
ba) states for antiferromagnetic wave vector Q= (vr,7r). 
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k^k+q 



^k4 



(6) 



with i,j = a, 5, El = -^(ei^) + (Aj^) , and Aj^ is the 
superconducting gap in the bonding (i — b) and anti- 
bonding (i — a) band. In the fohowing we assume that 
the pairing part of the Hamiltonian is symmetric with 
respect to the bilayers and given by 

Tipp = 5]A(k)(cl^^(k)cl^^(-k)+ct^^(k)4^(-k) 

k 

+h.c.) 

- E^W (cl^^(k)cl^,(-k) +ct^^(k)c[^(-k) 



+h.c.) , 



(7) 



where A(k) = -^(cosfca; — cos ky). It then follows that 
the pairing gap is the same for bonding and antibonding 
bands, implying A^ = A^ = Ak- However, the respec- 
tive Fermi surfaces in both bands are located at different 
momenta k in the Brillouin zone. In order to obtain the 
full x"^'" 1 we use the random-phase approximation (RPA). 
Within RPA the even and odd parts of the full spin sus- 
ceptibility are given by 



1 -5a(q)xo(q:'^) 



(8) 



where a = o,e and (7e,o(q) are the fermionic interaction 
vertices in the even and odd channels. 

We first consider the spin susceptibility at momenta 
close to Q = {t^,'^)- The dominant contribution to 
the susceptibilities comes from fermions near the hot 
spots, where both k and k + Q are close to the Fermi 
surface. In a dx2_y2--waYe superconductor with the 
above A(k) one has A(k-I-Q) = — A(k). As a con- 
sequence, Im^Q ° exhibits discontinuities due the open- 
ing of the superconducting gap2i. For the odd sus- 
ceptibility, Imxo'' and Imxo" exhibit a single disconti- 
nuity at ^fiQ) = |A^| + |A^^q|, where k is chosen 
such that £a(k) = eb(k + Q) = (see Fig[T]). Be- 
low this frequency, Imxg** = (at T = 0). At the 
same time, Imxg possesses two discontinuities located at 
nriQ) = |A£| + |A-+q| and Q'J'iQ) = |A^| + |A^+q|, 
where k is again chosen such that both fermions are at 
the Fermi surface (see Fig. [T]). Im^o is zero below the 
lower discontinuity, and jumps between two finite values 
at the higher discontinuity. Analyzing Eq. ([5]), we find 

Hence, in the even and odd channel, the susceptibility 
at low frequencies is purely real and, according to Eq.®, 
one finds that the bare Xo (Qi ^) (o^ — ^) behaves as 



(9) 



where Aq, > 0, and f{x) cx at small x, and f{x) ~ 
I log(l — a;)| near x = 1. Substituting this result in Eq.®, 
one finds that since f{x) changes between and oo when 
X changes between and 1, the susceptibilities in both 
the odd and even channel develop resonances below the 
thresholds of the particle-hole continuum, at frequencies 

fie.o where 1 = 5e,oXo °(Q: ^e,o)- 

As mentioned above, there are two reasons why the 
resonances in the even and odd channels occur at dif- 
ferent frequencies. One is that the even and odd free- 
fermion susceptibilities, Xo°{Qt^) are different, another 
one is that the interactions are different in the even and 
odd channels. The difference between the resonance fre- 
quencies in the even and odd channels due to the dif- 
ference in Xq°(Q,cj) arises predominantly from the fact 
that the (dimensionless) magnetic correlation lengths 
= 1 — ffcoXo °(Q: 0) are different in both channels al- 
ready in the normal state. Additional differences between 
Xo °(Q, 0) which arise in the superconducting state scale 
as Aq/ Ep, and are hence small and can be neglected. 
Assuming that the difference between the even and the 
odd resonance is small, and that the resonance frequen- 
cies are sufficiently low, such that f{x) in Eq.® scales 
as x'^, wc find after some simple algebra that 



no 



(10) 



The r.h.s. of the above equation is in turn related to the 
difference in the normal state static x via 



6-^ 



(11) 



The dominant contributions to the r.h.s. of Eq. (fTT|) come 
from fermions in hot regions near (0, tt) and (7r,0), for 
which the tj_ term in the dispersion, Eq.® reduces to 
±t±. Expanding the r.h.s. of Eq. (fTT|) to leading order in 
t±, we obtain 



^ Co ^ _ ,2 9^o 



= i1 



dujtPk 



(fife - iuj)'^{ek+q ~ iujf 



(12) 



where is the in-plane dispersion (i.e., Eq.® with t_\_ = 
0). Linearizing £k and Ek-i-Q in the hot regions as vp{kx + 
ky)/^/2 and vpikx — ky)/y/2^ respectively, substituting 
this expansion into the susceptibilities, and performing 
the integration, we obtain 



(13) 



where kmax ~ is the upper limit of the integration 
over momentum. Observe that the r.h.s. of Ea. p3p is 
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positive, implying that because of the difference in the 
even and odd susceptibiHties, the resonance in the even 
channel occurs at a larger frequency than the resonance 
in the odd channel. To estimate the strength of the ef- 
fect, we use vphp ^ O.beV ^ 2t, g = U/A, where U us 
the Hubbard interaction (in the RPA, TOJ7/27r ~ 1), set 
kmax = kp, and define J± = 4t\/U, J = W^/U. We 
then have 



47r J 



(14) 



The second source for the difference between fte and 
J7o is the difference in the interaction strength between 
the two channels. As mentioned above, within the t — J 
model, the two interactions are given by Jo,e = J\\ (q) ± 
J±. At q = Q, this effect alone leads to 



Jj_ 
J 



(15) 



where J = J|| (q = Q). We see that both effects described 
by Eqs.([T4]) and (llSp are in fact of the same order, and 
both lead to a larger fig compared to f^o- Moreover, the 
effect of the t± dependence of the interaction is larger, at 
least near optimal doping, where ~ 1- However, with 
decreasing doping, and hence increasing ^o, the role of 
the difference in the even and odd free-fermion suscepti- 
bilities may become more dominant. 

In Fig. [21 we present the results for the bare and 
full susceptibilities at optimal doping, corresponding to 
fi = — 1.195t, that were obtained from a numerical eval- 
uation of Eqs.® and ([8]). We see that Rexo in even and 




an observable difference between fie and flo- This agrees 
with our analytic treatment. Hence, the difference be- 
tween rig and arises from the difference in the effective 
interactions ge and go- To reproduce the experimentally 
measured frequency splitting between both resonances at 
Q, and the dispersion of the two modes (see below), we 
use 



5o.e(q) = Jo {1 - 0.1 [cos((?,) + cos(gy)]}±0.027Jo. (16) 

According to Eq.([T]), the first (second) term on the r.h.s. 
of the above equation can be interpreted as arising from 
the in-plane (out-of-plane) exchange interaction J|| (q) 
{J±)- 

We present the RPA susceptibilities Xrpa Q in 
Fig- [He). We see that both even and odd susceptibil- 
ities show resonance behavior. By construction, the res- 
onance in the even channel occurs at a larger frequency 
than the odd resonance. Accordingly, the intensity of 
the even resonance is smaller, which agrees well with the 
experimental observations^. 



; (a) odd mode 



(eV) 






1.4 



0-6 



FIG. 2: (color online) (a) Imxo'°, (b) ReXo'° and (c) Imx^p^ 
as a function of frequency at the antiferromagnetic wave vec- 
tor Q = 

odd channels are almost identical below 2Ao, i.e. the dif- 
ference in the susceptibilities is too small to give rise to 



FIG. 3: (color online) RPA results for magnetic excitations in 
a bilayered d^2_y2 superconductor at optimal doping. Calcu- 
lated Imx°(a) and Imx'^(b) obtained from Eq. ([8]) as a func- 
tion of momentum along q — rj{n,n)) and frequency in the 
SC state. 
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III. THE DISPERSION OF THE RESONANCE 
PEAK 



IV. DOPING DEPENDENCE OF THE EVEN 
AND ODD RESONANCES 



We next consider the dispersion of the even and 
the odd resonances and present in Fig. [3] a plot of 
Imx^P^ (q, f^) at optimal doping along q = r](Tr,Tr) as a 
function of frequency. The momentum dependence of the 
odd mode's frequency and intensity, shown in Fig.[3l^a), is 
quite similar to that of the resonance mode in the single- 
layer modeU^. In particular, away from Q three dis- 
continuities in Im^Q emerge, corresponding to scattering 
channels with momenta q, {2tt, 0) — q, and (2tt, 2tt) — q. 
The first momentum corresponds to a direct transition, 
while the last two momenta describe scattering processes 
involving Umklapp scatteringi^. As discussed before, the 
resonance can occur only at frequencies below the low- 
est discontinuity in Imx^^r^- Since the superconducting 
gap decreases towards the diagonal of the Brillouin Zone 
(BZ), the resonance dispersion follows the momentum de- 
pendence of the ph continuum, forming a parabolic-like 
shape^ii^. Upon reaching Qq ~ (0.8, 0.8)7r correspond- 
ing to the wave vector connecting the nodal points of the 
superconducting gap on the Fermi surface, the spin gap 
vanishes, and no resonance is possible. For even smaller 
q one finds that another resonance branch emerges, the 
so called Q* mode, arising from an umklapp transitio»i^. 

In contrast, the even part of the spin susceptibility ex- 
hibits six discontinuities in Im^o away from Q — (7r,7r). 
Intraband scattering within the bonding and antibond- 
ing bands each gives rise to three of these discontinuities. 
Similarly to the odd susceptibility, we find that a gen- 
uine resonance occurs only below the lowest discontinu- 
ity in Imxo' due to the direct transition with momentum 
q. This transition is again responsible for the parabolic- 
like shape of the even mode's dispersion, as shown in 
Fig. E^b). However, we find that the intensity of the 
even resonance falls off much faster as one moves away 
from Q than that of the odd one. Since the supercon- 
ducting gap and the splitting of the Fermi surfaces is 
zero along the diagonal of the BZ, the position of the 
so-called "silent band" is the same for the odd and for 
the even channel. Thus, both resonances merge together 
at Qo « (0.8,0.8)71- (see also FiglDJc)). Similar to the 
resonance in the odd channel, we find that for momenta 
smaller than Qq, an umklapp transition leads to the for- 
mation of a Q* mode in the even channel. However, its 
energy range is much smaller than that of the odd Q* 
mode due to the proximity to the ph continuum. For 
the odd as well as the even resonance mode we find that 
while the intensity of the Q mode (i.e., the mode origi- 
nating at Q) is largest along q = (tt, tjtt) and q ~ {rjir, tt), 
the Q* mode has its largest intensity along the diagonal 
direction, i.e., along q = 7y(7r,7r) and q = [(2 — r])^,!]^]. 
As previously discussecU^, this rotation of the intensity 
pattern by 45° reflects the qualitative difference in the 
origin of the first and second resonances. 



Next, we consider the doping dependence of the res- 
onance modes in the odd and even channels. In order 
to describe the doping dependence, it is necessary to 
know that of the superconducting gap as well as that 
of go.e(q). The doping dependence of the superconduct- 
ing gap, which is shown in Fig. IHb), is taken from recent 
ARPES experiments^ which suggest that the supercon- 
ducting gap increases by about 10 -f- 20% going from the 
optimally doped to the underdoped cuprates. In order to 
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FIG. 4: (color online) (a) Doping dependence of (a) the res- 
onance frequency at Q in the odd and even channel, and (b) 
the superconducting gap Ao and Jo. (c)-(f) Dispersion of the 
even and odd modes for various doping concentrations in the 
(c) underdoped, (d) optimally doped, and (e) and (f) over- 
doped regime. 

describe the doping dependence of (7o,e(q), we leave the 
momentum dependence of go,e(q) unchanged, and only 
change the overall prefactor, Jq, in Eq. (|16p . as a function 
of doping, by fitting the frequency of the resonance in 
the odd channel. The doping dependence of Jq is also 
shown in Fig. [IJb). We find that this procedure provides 
a satisfactory fit to the experimentally measured disper- 
sion of both resonance modes over a considerable range 
of doping. 

In Fig. [41(a) we present the doping dependence of the 
resonance in the even and odd channels at Q = (tt, tt). As 
expected from the discussions above, we find that with in- 
creasing doping, the energy splitting between both modes 
decreases, and for S = 0.21 is only about Aw^es ~ 1 
meV at Q, while for 5 = 0.15 one has Aw^es ~ 12 meV. 
This decrease in the splitting is observed over the entire 
dispersion of the resonance modes in the even and odd 
channels, which we present in Figs. |4l(c)-(f) for several 
different doping levels. In addition, we find that the dis- 
persion of the even mode exhibits a continuous downshift 
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with increasing doping, while that of the odd mode first 
shifts upwards with increasing doping in the underdoped 
systems, but shift downwards in the overdoped regime. 
In order to understand this quahtative difference between 
the underdoped and overdoped region, we note that in 
general, the doping dependence of the resonance modes 
is determined by that of the superconducting gap (which 
in turn determines that of the ph continuum) as well as 
that of 5^0,6 (q)- While a decrease of the superconducting 
gap, and hence a downward shift in frequency of the ph 
continuum leads to a downward shift of the resonances, a 
decrease of (?o,e(q), in contrast, leads to an upward shift 
of the modes' dispersion. 

Since the dispersion of the even resonance is located 
in frequency close to the ph continuum, and Rex3 varies 
strongly in the vicinity of the ph continuum due to its 
logarithmic singularity, it follows that the dispersion of 
the even resonance is rather insensitive to changes in 
(7e(q)- As a result, the doping dependence of the even 
resonance is predominantly determined by that of the 
ph continuum, exhibiting a continuous downward shift 
in energy with increasing doping. In contrast, in the un- 
derdoped regime, the energy difference between the ph 
continuum and the odd mode's dispersion is rather large, 
and Rcxq varies only weakly around the resonance fre- 
quency. As a result, the resonance frequency is very sen- 
sitive to changes in go{<l)- Therefore, it is the decrease 
in (7o(q) with increasing doping (and not the decrease 
in the superconducting gap) that determines the dop- 
ing dependence of the odd mode's dispersion and leads 
to its upward shift in energy in the underdoped regime. 
Around optimal doping, the odd mode's dispersion has 
become sufficiently close to the ph continuum, that the 
mode's further doping dependence is now determined by 
that of the ph continuum, and not any longer by that of 
(7o(q), similar to the case of the even mode. Hence, the 
two opposite effects arising from a decrease of the super- 
conducting gap and that of go(q) lead to the qualitatively 
different doping dependence of the odd mode's dispersion 
in the underdoped and overdoped regime. Note that with 
increasing doping, and the resulting downward shift of 
the ph continuum, the momentum range over which the 
Q*-mode can be observed, decreases. 

Finally, we briefly discuss the doping dependence of 
X'^'°{Q,uj ~ 0). If indeed, the suggested above, the odd 
and even resonance are transformed into the acoustic and 
optical branches of the spin wave dispersion in the anti- 
ferromagnetically ordered phase, one would expect that 
XoiQ,!-^ = 0) increases with decreasing doping. As a re- 
sult, one would see a downward shift in the odd mode's 
dispersion even for a doping independent Jq. One finds, 
however, that the doping dependence of Xq{QiO)j which 
is obtained from Eq.Q by simply changing the chemi- 
cal potential, /x, defies this expectation. This is shown 
in Fig[5l where we present the doping dependence of 
Xo'^(Q,0). Note that the even susceptibifity possesses 
two logarithmic divergences as a function of doping which 
occur when either the bonding or antibonding Fermi 




FIG. 5: (color online) ReXo'°(Q,0) as a function of doping 
concentration in the normal state. 



surfaces touches the van Hove (vH) points (±7r, 0) and 
(0, ±7r) and undergo a topological transition from a hole- 
like to and electron-like Fermi surface. These transitions 
occur at a doping level of a; sa 0.23 for the antibonding 
band and at x w 0.55 for the bonding band (not shown). 
In contrast, the odd susceptibility, which arises from scat- 
tering transitions between the bonding and antibonding 
bands does not exhibit a logarithmic divergence, but is 
simply enhanced and exhibits a finite maximum. If we 
define the minimum distance (in momentum space) of the 
bonding and antibonding Fermi surfaces to the vH point 
(0,7r) by kaifJ.) and ktiiJ,), respectively, then Rexo(Q,0) 
exhibits a maximum at that doping level for which the 
smaller of ka{n) and kb{p) possesses a maximum. Defin- 
ing fc,„m(/x) = Tam[ka{^.),kb{^)] one finds 



Rex°(Q, 0) ~ const. + ^— ^ arcsin 



(17) 



Note that for doping levels below that at which the 
van Hove singularity in Rexo(Q,0) or the maximmn in 
Rexo(Q, 0) occurs, the susceptibilities decrease monoton- 
ically with decreasing doping, as shown in Fig. [5l This 
doping dependence clearly reflects a shortcoming of the 
weak-coupling approach used above, which fails to cap- 
ture the strong correlation effects that are not only re- 
sponsible for the occurrence of antiferromagnetism, but 
are very likely also the key ingredients in the explanation 
of the pseudo-gap region in the underdoped cuprates. It 
is interesting to note in this context that recent studies of 
the doping dependence of xo(Q, 0) for a single layer sys- 
tem within the FLEX approach find that the vH singu- 
larity is eliminated by interaction effects, and that start- 
ing from the overdoped region Xo(Q5 0) increases mono- 
tonically with decreasing doping-22.. This shortcoming of 
the approach used above is effectively compensated by 
a phenomenologically introduced doping dependence of 
g'^'° which increases with decreasing doping. This phe- 
nomenological approach, however, does not allows us to 
fully explain the doping dependence of the resonant ex- 
citations in the underdoped cuprates. In particular, it 



7 



leaves open the question how the downward dispersion 
of the resonance mode observed in the optimally doped 
cuprates is transformed into the upward dispersion of the 
acoustic spin-wave branch. 

V. SUMMARY 

In this study, we have investigated the form of mag- 
netic resonance excitations in the even and odd spin 
channel of the bilayer cuprates in the superconducting 
state. We obtain a number of new results suggesting fur- 
ther experimental test that may finally resolve the long- 
standing question concerning the origin of the resonance 
peak. First, we show that the energy splitting between 
the even and odd resonances arises not only from a differ- 
ent interaction strength in both channels, but also from 
a the difference in the free-fermion susceptibilities in the 
even and odd channels. Both effects scale as ~ Jl/ J 
and lead to a frequency for the even resonance that is 
larger than that of the odd resonance. However, at least 
at optimal doping, the numerical prefactors are such that 
the energy splitting is dominated by the difference in the 
interaction strength and not by the difference in the free- 
fermion susceptibilities. Since the latter scales with 
the relative importance of these two effects might change 
in the underdoped cuprates. In agreement with previous 
resultsi^ii^ii^ we also find that the intensity of the even 
resonance is weaker than that of the odd resonance. Sec- 
ond, we computed the dispersion of the even resonance 
and showed that the even resonance also disperses down- 
ward as one moves away from Q = (tt, tt). Moreover, we 
demonstrated that the downward dispersion of the even 
mode is more parabolic than that of the odd channel. 
Third, we showed that there exists a second branch of 



the even resonance, similar to the recently observed sec- 
ond branch (the Q*-mode^- ) of the odd resonance, -i^^. 
We find, however, that in the even channel, this second 
branch is much narrower in energy than in the odd one. 
Fourth, we studied the doping dependence of the both 
resonance modes, and find that that of the even mode 
is determined by the downward shift of the ph contin- 
uum with increasing doping. In contrast, the upward 
shift in frequency of the odd resonance in the under- 
doped cuprates is determined by the decrease in Qo with 
increasing doping, while in the overdoped regime, the 
odd resonance follows the doping dependence of the ph 
continuum. Our results demonstrate that the structure 
of magnetic excitations in the superconducting state of 
the bilayered cuprates is dominated by the topology of 
the Fermi surface, the interaction strength in the even 
and odd channel, and the d^2 _y2-w<3Ne symmetry of the 
superconducting gap. 

Finally, we note that the experimental situation has 
recently been complicated by the report that an even 
resonance exists at incommensurate wave vectors only^. 
This result contradicts earlier studies which have found 
that the even resonance exhibits the largest intensity at 
Q = (tt, 7r)ii. The origin of this experimental discrep- 
ancy is currently unclear. 
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